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Abstract 

We explore gravitating relativistic spheres composed of an anisotropic, barotropic fluid. 
We assume a bi-polytropic equation of state which has a linear and a power-law terms. The 
generalized Tolman-Oppenheimer-Volkoff (TOV) equation which describes the hydrostatic 
equilibrium is obtained. The full system of equations are solved for solutions which are 
regular at the origin and asymptotically flat. Conditions for the appearance of horizon and a 
basic treatment of stability are also discussed. 


1 Introduction 

Though very pioneering, the Lane-Emden equation is still a versatile and useful modeling device 
in astrophysics [T]. To obtain the Lane-Emden equation one needs to invoke three essential 
basics: spherical symmetry, hydrostatics equilibrium, and a polytropic equation of state [2]. Of 
course, the relativistic version of the Lane-Emden equation is now at our disposal [3]. In general 
relativity one needs to replace the classical hydrostatic equilibrium equation with its relativistic 
dual. The well-known TOV (Tolman-Oppenheimer-Volkoff) equation is in fact the relativistic 
version of hydrostatic equilibrium HElIlllZ!. Generally, the Lane-Emden equation can be derived 
using the polytropic equation of state in the form P = with n being called the polytropic 

index. The rationale for postulating such a relation between density and pressure originates 
from thermodynamics. In thermodynamics, a polytropic equation comes from the adiabatic 
condition dQ = 0 where dQ is the transferred heat. One of the immediate consequences then is 
Pp^ = constant in which P is the pressure, p is the density and 7 = the ratio of the two 
specific heats [8]. The isothermal and the adiabatic processes are the most important examples 
of a polytropic process. To find a comprehensive analysis about the non-relativistic version of 
the Lane-Emden equation and its solutions one can refer to Clayton (1968) [9]. 

Assuming a polytropic star makes the relativistic structure analysis easier and more inter¬ 
pretable. Solutions should satisfy appropriate boundary condition(s) in order to explain various 
astrophysical objects from compact stars to galactic objects [TO]. Almost all polytropic models 
assume an isotropic pressure. Here, we allow anisotropic cases by assuming different tangen¬ 
tial and radial pressure components in the energy-momentum tensor. All principal pressures 
(stresses) are assumed to have double-polytropic equations of state [TT|. A number of authors 
have investigated anisotropic models (see, [I2]-[I7]). Assuming anisotropic stresses is by no means 
un-natural. For example, in a compact star, although the radial pressure vanishes at the surface, 
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one still could postulate a tangential pressure to exist. While the latter does not alter the spher¬ 
ical symmetry, it may create some streaming fluid motions. The Reissner-Nordstrom solution 
is another relativistic solution which is supported by anisotropic stresses. A charged black hole 
readily results in negative radial pressure, while the tangential pressure remains positive. On 
galactic scales, a spherical galaxy is well described by the King-Michie family of models [isms], 
which involves an anisotropic pressure. This anisotropic pressure originates from the fact the 
phase space distribution function of the stars depends on coordinates and velocities only through 
the constants of motion (energy and angular momentum). If the distribution function depends 
on the magnitude of the angular momentum, this will lead to an anisotropic pressure, even if the 
system is spherically symmetric. 

In what follows, we first derive the TOV equation for an anisotropic fluid with spherical 
symmetry (next section). In section [3l we will derive analytical solutions with regular behavior 
at the center. Section 0] is devoted to a study of Killing horizons and energy conditions. The last 
section reflects concluding remarks. 


2 TOV equation 

Although for most cases the classical hydrostatic equation suffices, when the gravitational field 
becomes strong, using a general relativistic version of the hydrostatic equilibrium equation be¬ 
comes inevitable. If one starts from a spherically symmetric metric assuming an isotropic energy- 
momentum tensor one ends with the famous TOV equation [20] 

dP _ Gjpc^ + P){M{r) + Anr^P) 
dr r‘^{c^ — 2GM{r)/r) 

where P is the pressure and M(r) is the mass variable defined by 

pr 

M{r) = / p(/)47r/^dr'. (2) 

Jo 

From now on, we release the assumption of isotropic pressure and let the tangential and radial 
pressures to be different. Let us take the following static and spherically symmetric metric for 
describing a gravitating relativistic sphere in Schwarzchild coordinates: = {t, r, 6, (p) 

ds^ = (—1 -|- /(r))c^dt^ + (1 “ /(^))~^dr^ -|- r^dO^, (3) 


where /(r) is a function of the radial coordinate r. Introducing an orthonormal spherically 
symmetric basis 

9^ = {—1 + f{r))^^‘^cdt , 0^ = (1 —/(r))“^/^dr, 6‘^ = rd'd, 6^ = r sin i!) d(f>, (4) 

and using dO^ = —v for the torsion-free spacetime, the nonzero components of the 
connection form will be 


0 _ 1 _ fi’’’) flO 

to I — U) 0 — 


2(l-/(r))V2 ’ 


,.2 ,.i ( l -/( 0)^/%2 

a; 1 — —u) 2 — -V ) 


,3 _ _(l-/(0)^/%3 


OJ 1 = -UJ 3 = 




(5) 

( 6 ) 
(7) 
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( 8 ) 


.3 .2 COt??„3 


U} 2 = -W 3 = 


e\ 


We also have the energy-momentum conservation law in the form 

D*To, = 0, (9) 

where D stands for absolute exterior differential m and Ta is the energy-momentum 1-form 

Ta = ( 10 ) 

In our case, *Ta can be written in the form 

( 11 ) 

Note that there is no summation over indices in (|lip . and q°‘ and are defined as 

= diag(-pc^ Pr, Pt, Pt), 

and 


r?" = 

where Pr and Pt are the radial and tangential pressures (stresses). 

From the conservation law Q one obtains 

dq°‘ A ry" = — ^ A r]^{q^ — q°‘), 

0 

for a = 1 we have 

dPr A A r]^{q^ - q^) + oj^ 2 - g^) + wS A rf’{q^ - q^). 

Substituting for the connection forms from Q-dS]) we obtain 

(1 - /(r))V2p'(r)0i A ^ 


( 12 ) 

(13) 


(14) 


(15) 


(1 - /(r))i/2^_2 


{pc^ — Pt{r)){9‘^ Ap'^ +9^ A p^). 


(16) 


Using the metric Q and matter field ()12p . the first two components of the Einstein field equations 
reduce to 

1 1 

(17) 


1 


GOq = G\ = --fir) - A, fir) 


which is equivalent to Pr = —pc?. On the other hand (I17p can be written as 


/ «/ SttG 2 
r/ r ) = -^pr\ 

where prime stands for differentiation with respect to r. Solving for /(r) yields 


f{r) = - + [ AtcG pr'"^drf 

r rc^ J 


(18) 


(19) 
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with b is the constant of integration. Though having a singularity (black hole) at the origin is 
among the possibilities and should be explored separately, for the moment we ignore the first 
term of the above equation in the following relations. Consequently /(r) reduces to 


f{r) = 


2GM{r) 

c^r 


where M(r) is the gravitational mass inside radius r in such a way that 

pr 

M{r) = / Anr'^pdr'. 

Jo 

Substituting for /(r) and f'{r) in (fTUD we hnally obtain 


dPr _ G{pc^ + Pr){M{r)c^ + Ai^r^Pr) 2{Pr - Pt) 
dr r‘^{c^ — 2GM{r)/r) r 


( 20 ) 


( 21 ) 


( 22 ) 


which is the generalized Tolman-Oppenheimer-Volkoff (TOV) equation. It can be seen that Eq. 
()22p readily reduces to the standard TOV ([T]) when P^ = Pt- We will look for solutions of the 
generalized TOV equation which are asymptotically Schwarzschild, in the following section. 


3 Analytical solutions 

In this section, we will derive analytical solutions of the generalized TOV equation which are 
non-singular throughout the space and asymptotically approach the Schwarzschild spacetime. In 
what follows, we assume that the radial and tangential pressures obey the following barotropic 
equations of state (EoS): 

n 

Pr = wpc^ + (23) 

Pc 

m 

Pt = Wipc^ + (24) 

Pc 

here pc is the central density and w,w',wi,W2 are the dimensionless EoS parameters. These 
EoSs are composed of two terms. The hrst term is linear which is the most widely used one in 
astrophysics and cosmology. The second term allows a polytropic behavior which is familiar in 
astrophysical applications, leading to the Lane-Emden equation in the theory of stellar structure 
[9]. Since from Eq. (|17p we have the equality Pr = —pc^ one obtains w = —1 and u)' = 0 in Eq. 
(I23p . This value for w shows that we are encountering a dark-energy-like situation. Existence 
of negative pressure warns us about possible violation of energy conditions. However, nowadays 
energy condition violating matter has become much more popular in the literature (e.g. dark 
energy, wormhole theory, quantum effects). In particular, the existence of a strongly negative 
pressure is an essential ingredient in dark energy models. Lobo [22] presented a generalization 
of the gravastar picture, by considering a matching of an interior solution governed by the dark 
energy equation of state re < 1/3, to an exterior Schwarzschild vacuum solution at a junction 
interface. He assumed an isotropic pressure there. The present work differs from [22] in three 
respects. First, the fluid pressure is assumed to be anisotropic. Second, the EoS is not assumed 
to be linear, and third, we do not cut off the matter density at a certain radius, thus there is no 
need to apply matching conditions at the surface. 

The 22 and 33 components of the Einstein equations give 

G\ = G\ = -^r{r)-^f{r). (25) 
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( 26 ) 


Eqs. CZI) and ()25l) can be simplified by introducing u = ^ in the following form 

3 SttG 

- u-^fu + u^f = 


r .2 


4 ^ 

u tuu= - 


IGttG 


-Pt. 


(27) 


respectively. Here, the subscript u indicates derivative with respect to the variable u. Taking 
derivative with respect to u from Eq. (j26p and then multiplying by u and finally using Eq. (I2ZD, 
one obtains 

p^{u) 1 dp{u) 


p{u){l + Wl) + W2 


which can be solved for density p{r): 

p{r) = 


2dlnu'' 


Pc 


JLp(l+«)i)(m-l) _ Q,]m-i ’ 


where 


f3= {xq ™ + a), xo = 


pjro) 

Pc 


a = 


W2 


1 


II , ro = —. 
1 + rci Mo 


(28) 


(29) 


(30) 


Eq. (|29p is the general form of the density in all of the models which will be discussed separately 
in the following sub-sections. At the center of the star r = 0 the boundary condition p{r) = pc 
implies = 1 corresponding to a = ±1 if 1 -|- tci > 0. 

3.1 Case 1 


The first case corresponds to m = - in Eq. (|2^ . By substituting the values a = —1 and W 2 = —2 

into Eq. (1301) . we obtain wi = 1. Now the tangential pressure of the anisotropic star for this 
special case becomes 


P, = PC^ - 2 

\/1p 

the mass function can be obtained by performing the integration in Eq. m 

r 


M{r) = 


2-KpcrQ- 




arctan(V^(-)) - 
Mo 


Mo 


2/3(/3$ + 1) 


and the function /(r) becomes 


/(m) = 


2GM{r) AnGpcro" 


c^r 




arctan(V^(^)) - 
Mo 


r 

ro 


2/3(/3^ + l) 


The total gravitational mass of the system will be 

M = lim M(r) = 


/3i ^ 


(31) 


(32) 


(33) 


(34) 


therefore, the metric is asymptotically flat (/ —>• 0). Furthermore, at the center the mass func¬ 
tion becomes zero, as expected for a non-singular model. Thus this case satisfies the physical 
requirements of a non-singular relativistic gravitating system. 

The metric and density functions are shown in Figures [T] and [21 respectively. 
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Figure 1: The metric function goo = f{x) — 1 as a function of x = r/ro is plotted for Case I. Note that 
there are two coordinate singularities (Killing horizons) for the values C = 2 and /3 = 0.5. For (7=1 
and /3 = 0.7, there are no Killing horizons. For (7 = 1 and /3 = 0.5 there are no Killing horizons a naked 
singularity exists at x = 0. 



r 

Figure 2: Energy density p{r) as a function of radial distance r is depicted for case 1. 


3.2 Case 2 

Setting a = — 1 and W 2 = —3 we get wi = 2, the tangential equation of state becomes 


Pt = 2pc^ - 3 ^, 

by substituting W 2 in Eq. (p0]l . the mass function is obtained as 

3 (JLp 


M(r) = 


3 l + 


furthermore, 


_ 2GM{r) _ SnGpcro^ 
•^7 j — ^2 


c^r 3c^r 1 + / 3 (^)^ 


(35) 


(36) 


(37) 
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The total gravitational mass of the sphere is 


M = lim M(r) = . (38) 

r—>oo 3p 

The metric for this case is asymptotically flat and at the center of star the mass function equals 
zero. The metric and density functions for this case are shown in Figures [3] and 01 respectively. 



Figure 3: The metric function goo as a function of the dimensionless radius for case 2. For (7=1 
and P = 0.2, there are two Killing horizons. For (7 = 1 and /3 = 0.38 the two Killing horizons become 
degenerate. Finally, with (7 = 1 and /3 = 0.5 there are no Killing horizons. This case does not have any 
curvature singularities. 



r 

Figure 4: Energy density p(r) as a function of radial distance r is plotted for case 2. 


3.3 Case 3 

Setting a = —1 and W 2 = —6 we have wi = 5, and the tangential equation of state becomes 

Pt = 5pc^-6^, (39) 
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the mass function is obtained as 


M(r) = 


TT/Ocr-O" 


ro 


1+Pi-r 

ro 


+ 


arctan y/p (—)’ 

r^ 


and 


, 2GM{r) 2TTGpcro^ 
/W = —3— -- 


c^r 


3c^7 


^arctan v^( — )^ 


(-) 


The total gravitational mass is obtained as 


M = 


l + /3(-)6 

ro 


TT^Pcrp^ 


+ 


ro 


VP 


(40) 


(41) 


(42) 


Therefore, the metric is asymptotically flat and at the center of the system the mass function 
vanishes. The metric and density functions for this case are shown in Figures [5] and [6l respectively. 



Figure 5: Two Killing horizons and a geometric singularity at a: = 0 for case 3, corresponding to C = 2 
and /3 = 0.1. C = 1.45 and (3 = 0.1 leads to degenerate Killing horizons. For C = l,/3 = 1-5, there are no 
Killing horizons but a naked singularity exists at a; = 0. 



r 


Figure 6: Energy density p(r) as a function of radial distance r for case 3. 
























3.4 Case 4 


When a = —1 and W 2 = —1 then lei = 0 and the tangential pressures becomes 


R = - 


p2C^ 

7^’ 


and the mass function will be given by 

47rpcro^ 


M(r) = 


/?2 


21n(/3- + l) 
r 


1 


ro 




+ 


1 


/3(/3- + 1) 
ro 


and 


2GM{r) S-nGpcVo^ 

RR = —:2— =- 


c^r 


13^ Gr 


21n(/3—+ 1) 
r 


1 


ro 


/? 


+ 


1 


/?(/?- + 1) /? 


(43) 


(44) 


(45) 


However, the total mass diverges in this case. In order to have physically viable solution, the 
total mass should be finite. Moreover, we have a singular mass at the origin. 


3.5 Case 5 

3 1 

When a = — 1 and W 2 = — then wi = — and the tangential pressure can be expressed as 




Pt = 


pc 3/9 2 C' 


“^y/Pc 


(46) 


The mass function will be 


... . STTpcro^ 

= -w~ 


- f —5- + ln(l + / 3 (—)2) — 1 

! + /?(-)§ ^0 


ro 


(47) 


and 


2GM{r) IQi^Gpcro^ 

/(O = . 2 . “- 


c^r 


3^2 


l + /3(-)^ 
ro 


+ ln(l + /3(—)i)-l 

ro 


In this case, too, we do not have a regular behavior at the infinity and the center, since 

r —)■ 00 , M —00 and / —)■ 0, 

IQiipcro^ 


0, M 


3/32 


and / —> 00 . 


3 5 

We have also reached analytical solutions for (rci = -,W 2 = —-) and (rci = 3,W2 
However, the solutions are too lengthy to be shown here. 


(48) 

(49) 

(50) 
-4). 


9 



4 Killing Horizons and Energy Conditions 


Since the spacetimes considered here are all static, we have the Killing vector K = dt and Killing 
horizon for the above cases can be obtained by putting the time-time component of the metric ([3]) 
equal to zero which gives /{xh) = 1, where xh = Unfortunately, analytical solutions for rjj 
could not be found. Therefore, we demonstrate qualitative f{x) — 1 behavior of the function for 
different values of C and {3 by plotting the cases mentioned in the previous section (see Figures 

([I])-®. Note that we have defined C = for convenience. 

As it is seen in Fig. ([T|) which is plotted for Case 1, the function f(x) — 1 has two killing 
horizons for the values (7 = 2 and (3 = 0.5 with a singularity at x = 0. With (7 = 1 and (3 = 0.7 
there is no Killing horizon. For (7 = 1 and (3 = 0.5 there are no Killing horizons but there is 
a singularity at x = 0. Fig. ([3]) is plotted for Case 2. With (7 = 1 and /3 = 0.2, there are two 
Killing horizons. For (7 = 1 and /3 = 0.38 a degenerate root exists which is regarded as the 
extremal Killing horizon. Finally, with (7 = 1 and /3 = 0.5 there is no Killing horizon. It should 
be noted that this case does not have any geometrical singularity. 

In Fig. ([5]), there are two Killing horizons with a singularity at x = 0, corresponding to 
(7 = 2 and (3 = 0.1 are shown. (7 = 1.45 and (3 = 0.1 leads to an extremal Killing horizon. For 
(7 = 1, /3 = 1.5 we do not have any Killing horizon but there is a singularity at x = 0. 

The metric function is not plotted for cases 4 and 5, since they are not physically interesting. 

Here, we examine weak energy condition (WEC) and strong energy conditions (SEC) for 
the cases introduced in the previous section. WEC requires T^i.V^V'^ > 0 for any non-spacelike 
vector field which leads to |18j 


P > 0, 

1 , 


P + 77 . > 0, p + Tt ^ 0. 


(51) 


SEC states that TnyV^^V^ > -T\V"Va for any timelike vector U^, or 


p + 77 . > 0, p + T) > 0, p -\- Pr ~\- ‘iPt P 0. 


(52) 


For Case 1, Eq. (I3TD reduces to p > 0 and p < Pc which are both satisfied since the density is 
maximum at the center and the density is assumed to be positive everywhere. Therefore, WEC 
is fully satisfied. SEC leads to pc > 4p which is also satisfied. Thus this case satisfies the physical 
requirements of a realistic relativistic model. 

WEC is satisfied for Case 2, similar to the previous case. SEC is also satisfied, provided that 
9 

Pc > -p. This case also satisfies the physical properties that should be satisfied by a realistic 

model. For the remaining cases, the WEC remains satisfied as before. SEC for Case 3, implies 
36 

Pc > —p. However, for Case 4 we obtain p < 0 which contradicts positive energy theorem. 
25 

Finally, for the last case, SEC requires pc > 9p. 

5 Stability and Horizon 

As it was shown, there appear Killing Horizons for some cases. The surface gravity k associated 
with the time-like Killing vector K = dt can be obtained via m 


1 


= -- {dK, dK) \h. 


(53) 
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K can be calculated for all of the cases considered before. As an example, for the second case we 
have 




2rH 


3/3r; 


l + /3(—)3 r3(l+ /3(—)3)2 
ro ro 


(54) 


The area of the outer Killing horizon can be found straightforwardly by setting r = rjj and 
dt = dr = 0 in Eq. ([3|) yielding 

ds^ = r|^dn^. (55) 

The horizon area is then the integral of the induced volume element, 

^ = / \/M dd d(/>, (56) 


so the area of the Killing horizon is simply 


A = ATrrjj. (57) 

Therefore, the area of the Killing horizon for all cases is given by the usual black hole area. The 
question whether the derived solutions are stable or not requires a separate full study. However, 
as a first check, we have plotted the total mass versus central density for the case 1 to case 3 
models. It can be seen from these figures that the total gravitational mass has a minimum at 
certain central densities, signalling gross stability. 



Figure 7: Total mass as a function of central density for case 1 model (all parameters are set to one). 
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Figure 8: Total mass as a function of central density for case 2 model (all parameters are set to one). 



Figure 9: Total mass as a function of central density for case 3 model (all parameters are set to one). 
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Note that when there are two horizons, the coordinates t and r become spacelike and timelike, 
respectively, in the region between the two horizons. Thus the situation is somehow similar to 
the non-extremal Reissner-Nordstrom black hole. 

6 Summary and Conclusion 

In this paper, we started with generalizing the TOV equation for a gravitating relativistic sphere 
with an anisotropic, barotropic fluid. We adapted an equation of state which has a linear term 
plus a power-law term which is encountered in various polytropic fluids. The full system of equa¬ 
tions was then solved exactly for solutions which are regular at r = 0 and smoothly approach 
an asymptotically flat spacetime. Some singular and asymptotically non-flat cases were also 
mentioned, without going into details. Finally, conditions for the emergence of horizons were 
examined. It was shown that in most cases, it is possible to avoid horizons, have two horizons, 
and also have a single degenerate horizon. The status of weak and strong energy conditions were 
also discussed. Finally, case 2 and case 3 models were shown to be stable, in the sense that the 
total gravitational mass is minimized for certain central densities. 
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